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Based on the Lippmann-Schwinger equation approach, a generalized Liischer's formula in 1 + 1 
dimensions for two particles scattering in both the elastic and coupled-channel cases in moving 
frames is derived. A 2D coupled-channel scattering lattice model is presented, the model represents 
a two-coupled-channel resonant scattering scalars system. The Monte Carlo simulation is performed 
on finite lattices and in various moving frames. The 2D generalized Liischer's fornrula is used to 
extract the scattering amplitudes for the coupled-channel system from the discrete finite-volume 
spectrum. 
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I. INTRODUCTION 

In recent years, remarkable progresses have been made 
on hadrons scattering in lattice QCD from both the the- 
oretical algorithm of extracting scattering amplitudes 
from lattice data [Tl-fT2] and the practical lattice QCD 
computational algorithm aspect [T51 - H7] . Since Luscher 
proposed the elastic scattering formalism in a finite vol- 
ume pQ, the framework has been quickly extended to 
moving frames [2H5], and to coupled-channel scatter- 
ing [7l412|. The finite volume scattering formalism has 
been successfully used by the lattice community to ex- 
tract elastic hadron-hadron scattering phase shifts [T51 - 
|2B]. Realistic lattice QCD computations on coupled- 
channel hadron-hadron scattering are under way. 

For the purpose of demonstrating the feasibility of ex- 
tracting coupled-channel scattering amplitudes from lat- 
tice data and discussing some issues, such as, finite size 
effects, in this work, we present a coupled-channel scat- 
tering lattice model in 2D. Our model is a direct general- 
ization of a 2D single channel scattering lattice model in 
|2"T] . The advantage of scattering in 2D is that only finite 
numbers of scattering amplitudes contribute in one spa- 
tial dimensional scattering theory, and the relation be- 
tween phase shift and energy level in Liischer's formula 
in 2D [T3] appears more transparent. Our 2D lattice 
model represents a coupled-channel resonant scattering 
system with three species of scalar fields (</>, a, p), where 
the scalar field p acts as a resonance which couples to 
both 2(j> and 2a channels. The Monte Carlo simulation 
are carried out on various lattice sizes and in different 
moving frames. We also present the derivation of 2D 
Liischer's formulae in a general moving frame and for 
a coupled-channel system. The derivation is based on 
the Lippmann-Schwinger equation approach presented in 
[12] . These formulae are used in the end to extract the 
scattering amplitudes from Monte Carlo simulation data. 
The finite size effect on extracting scattering amplitudes 



(phase shifts and inelasticity) from lattice data is also 
addressed in this work. Although, our model is formu- 
lated and computed in 1 + 1 dimensions, it still captures 
many of the features of hadrons scattering in a real 3 + 1 
dimensional QCD computation, and sheds some light on 
the future coupled-channel hadron-hadron scattering lat- 
tice QCD calculation. 

The paper is organized as follows. A discussion of elas- 
tic scattering in a finite volume is given in Sectionln] with 



extension to the coupled channel system in Section III 
The 2D lattice model, the Monte Carlo simulation and 
data analysis are described in Section |IV| The summary 
and outlook are given in Section [Vl 



II. LUSCHER'S FORMULA IN 1 + 1 
DIMENSIONS 

For completeness, we first present the basic scattering 
theory in 1 + 1 dimensions. Based on the Lippmann- 
Schwinger equation approach, a generalized Liischer's 
formula in 1 + 1 dimensions for two particles elastic scat- 
tering in moving franres is presented in the end of this 
Section. 



A. Two-particle scattering in infinite volume 

We consider spinless particles scattering in a symmet- 
ric potential V(—x) — V(x), the mass of scalar parti- 
cles is m. The wave function of scattering particles in 
center of mass frame satisfies the relativistic Lippmann- 
Schwinger equation 



ip(x) 



dx Gq{x — x; i/s)V(x')ip(x'), (1) 



where the center of mass frame energy is yfs and the 
free-particle Green's function is given by 
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G Q (x;y/s) 



dq 



,«/.!■ 



-oo 2tt ^s - 2 v / q 2 + m 2 ' 



(2) 



The Green's function can be further written as a oscil- 
lating term and an exponentially decaying term over the 
separation of two particles. The singularities of inte- 
grand in Eq.pl) on the complex q plane are two poles on 
real axis q — ±fc and two branch cuts on imaginary axis 
± [im, ioo], see Figfl] Therefore, for x > 0, we choose the 
contour C\ + C'2 to include pole q = k and cross the cut 
[im,ioo], and for x < 0, we choose the contour C\ + C3 
to include pole q — —k and cross the cut — [ira, ioo], as 
shown in Figfl] Thus, contour integral leads to 



G (x, y/s) 



-%- — e 
4fc 



ik\x\ 



2tt 



-p\x\ 



k 2 



(3) 

where k = vs-4m j g momentum of particle in CM frame. 
At large separations, the free Green's function can be 
approximated by the oscillating term only 



G (x-x',\/s) 

\A>\*'\ .\/s 



__ ik\x\ y^ 

4fc ^ 

v=± 



Yp(x)Y v (x r )J v (kx'), (4) 



where the functions Yp{x) and J-p(kx) are defined by 



Y + (x) = 1, Y_(x) 



x 

R 



(5) 



J+{kx) = cos k\x\, J^(kx) = isink\x\. (6) 

Such that Yp(x) and J-p(kx) resemble the spherical har- 
monic and Bessel functions in three spatial dimensions, 
and Yp(x) is the parity eigenstate with eigenvalue V. 
The continuous rotation symmetry in three dimensions 
reduces to discrete spatial reflection x — > — x in one 
spatial dimension, thus, the partial wave expansion of 
wave function in three dimensions reduce to the expan- 
sion of the wave function in terms of parity eigenstates 
i'ix) — J2v=± cvi^vix), where ip-p(-x) — Vtp-pix). 

For a potential V which falls at large separations, 
Eq. (IT]) is solved outside the range of the potential by 



ip(x) 



\x\>R 



J2 c v Y v {x) \j v (kx) + ie tklxl /p(fc)l , (7) 



v=± 



where R denotes to the effective range of potential and 
the free solutions has been also included in Eq.([7|. The 
scattering amplitudes are defined by 



cvfv(k) 



4fc 



dx'Yp(x')J£(kx')V(x')ilj{x'), (8) 



which up to the inelastic threshold can be parametrized 
by scattering phase shift 



f v (k) =e idv sin5 v 



(9) 




FIG. 1: The integration contours and singularities of free 
Green's function in Eq.S on complex q plane. 



B. Two-particle scattering on a torus 

Now we consider the theory in a one spatial dimen- 
sional box with periodic boundary conditions. In lattice 
QCD calculations, the computations are usually done in 
the moving frame of the two-particle system [2]. After 
the system is boosted back to the CM frame, the shape 
of cubic box in moving frame is deformed in CM frame 
due to Lorentz contraction. Similarly, in the one spatial 
dimension, the volume of a one dimensional box, L, in 
a moving frame with total momentum P = ^-d,d G Z 

becomes jL in CM frame, where 7 = Jl + ^- is the 
Lorentz contraction factor. 

Taking into account the Lorentz contraction effect as 
well, we divide the integral over x' into a sum of integrals 
over each translated box in Eq.pJ), giving, 



^ L \x) 



E 



dx'Go(x — x' — 77J-L; y/s) 



x V(x' + jnL)^ L) (x' + jnL). 



(10) 



The wave function in CM frame satisfies the boundary 
condition |2] of 



^ L \x + jnL) = e l ^ nL ^ L \x), 



(11) 



Using the periodicity of the potential V(x' + 771L) 
V(x'), we have 



^ p \x) 



-_A 



dx'G P (x - x'; y/s)V(x')^ L ' p \x'), 

where the periodic Green's function is given by 
G P (x - x': Vs) = J2 G °( x ~ x '~ T" L; \^)e^ nL . (13) 



By using the Poisson summation form ula 



hY. n& ^ nL = ^£„ e z<*(£ + f£n), Eq.(|l3|) 



2tt ^m£Z ^ 7L 

can be reexpressed as 

G p {x-x i ;V~s) = \Y 



iq{x—x') 



7^iX V^-2Vg 2 + m 2, 



(14) 



where P d = {g e M|<7 = ff(n + f),for u6Z}. 

As in the infinite volume case, the periodic Green's 
function Eq. ( 14 ) can be shown to consist of an oscillatory 



part and an exponentially decaying part which can be 
neglected for large volume 7L > m~ l . The remaining 
oscillatory part takes the form 



G P {x ~ x'; v/i) -> -i^- Yl e ik \ x -*-T nL \e i $ nL , 

nez 



(15) 



where we have used Eq.p| and Eq.p"3"|). The infinite 
sum in Eq.(15) can be done analytically, the details are 
presented in Appendix |A"{ so that 



Gp(x-x';Vs) 

M>|z'l .\fs 
~ ~*4fc 



Y Y v (x)Y v (x')J v (kx') 



v=± 



^N-("l-,cot 7fc£ 2 +7rd 



J-p(kx) 



(16) 



Using the definition of scattering amplitudes in Eq.Q, 
we can express the wave function as 



jkL + nd 



v=± 

k\x\ 



ic 



i + cot 



J-p(kx) 



(17) 



Matching the wave function in finite box given by 
Eq.(17l to the wave function in infinite volume given by 
Eq.(7Tat a arbitrary \x\ > R, we obtain 



Y c vY v (x)Jv(kx)f v (k) 
V=± 



1 7/cL + nd 

+ i + cot 



Jv{k) ■ • " 2 

which has non-trivial solution when 

jkL + nd 
cot o-p + cot = 0. 



(18) 



(19) 



III. COUPLED-CHANNEL SCATTERING IN 
1 + 1 DIMENSIONS 



The previous discussion in sectionlTTlcan be generalized 
to a coupled channel system by including another species 
of scalar fields, let's name two species of particles <fi and 
a, the masses are m$, a . The coupled channel wave func- 
tion has the form of i(>(x) = ^2 a=(t> a ip a (x), and i/j a (x) 
satisfies equation 



/oo 
dx'G^x-x'-JS) Y V a p(x')^(x'), 

■°° P=<t>,a 

. ( 20 ) 

A 2 x 2 matrix of coupled-channel scattering amplitudes 
can be defined by 



^>fv=-J& / dx'Y v (x')J r (k a x')V a0 (x')^(x'), 

4K Q J-00 

(21) 

where k a — — — ^ — - is the CM frame scattering mo- 
mentum in channel a. Neglecting exponentially decaying 
terms and also include the free solution, we have for the 
wave function in channel a 



r(x 



\x\>R 



£*w 



v=± 



c^Mk a x) + ie ik ^Y4fr 



(22) 



Extending the single channel derivation in finite- 
volume to the two-channel system, one obtains, 



P=± 13 



(23) 



Matching the wave function in finite- volume, Eq.(23l to 



the wave function in infinite volume, Eq.(22), we can 
derive a condition for non-trivial solutions 



1 ik^L + ird 

—rT + i + cot — ^—- 

r<t><t> 9, 




1 lk„L + nd 
+ i + cot 



cot 



jk^L + ird 



cot 



iKL 



ird\ (fv) 



ft>9 fa 
J-p Jp 



lp Jv 



(24) 



The scattering amplitudes can be parametrized by three 
real parameters: two phase shifts 5%, and an inelasticity 
Vv, 



}v ~ 2i 



1 



p<\8 



rap _ 

Jv — 



yr^e^+o 



1 Jv 



Thus, we can also write the Eq.(24) as 



cos ( 7-L 



KA, -\-K a 



cos ( 7^^ — 



sf, + s v 



6*,-6 v 



(25) 



(26) 



IV. 



THE ISING MODEL FOR COUPLED 
CHANNEL SCATTERING 



To simulate a coupled channel scattering system, we 
build a model with two light mass fields (0, a) coupled 
to a heavier mass field p with two 3-point couplings, pcj) 2 
and pa 2 . The physical masses of the fields are calibrated 
to be at the region 2m^ < 2m a < m p < Am^. For elastic 
scattering, the Ising model has been used and tested in 
both 1 + 1 [57] and 3 + 1 [5] dimensions by coupling two 
Ising fields, (f> and p, together through a 3-point nonlocal 
interaction. For our purpose, we could introduce one 
more species of Ising field a and another 3-point term to 
couple a and p together, where p field gives rise to the 
resonant behavior in both (fxf> and aa channels. 

The action is given by 

S = — 2_. K a / j ce{x)a{x + p.) 

a—<p.a,p x.p 

+ 2 ffpw5>(^)/3(20/^ + A), (27) 

where x — (xq,x\) is coordinates of Euclidean T x L 
lattice site and ft denotes the unit vector in direction 
p. The values of the fields are restricted to ±1, and the 
periodic boundary condition has been applied in Monte 
Carlo simulation. In the scaling limit, the Ising model 
represent a lattice (j) A theory, thus, the action in Eq.(27) 
effectively describes an interacting theory of [2] 






d A x 



cl 2 j 



- {da) + -m a a + —a 



9 pa a 2 



(28) 



in Euclidean space. By adjusting the masses and cou- 
pling constants, we could have an resonance p sit above 



both 2cf> and 2a thresholds, and couple to both channels 
by interaction terms g p ^p4> 2 and gpaapo 2 respectively. 
Therefore, the lattice Monte Carlo simulation by using 
the action in Eq.(27l is expected to imitate a coupled- 
channel scattering model: (jxp + aa HpH cf>(f> + aa. Due 
to the Bose-symmetry, only scattering amplitudes with 
positive parity contribute in this model. 



Cluster algorithm for coupled-channel Ising 
model 



An generalized cluster algorithm is used in our simu- 
lation, similar to the cluster algorithm in [27J , we update 
p, <f> and a fields alternately. 

Updating the p field: Bonds between neighbored spins 
of equal sign are kept with the probability 1 — e~ 2Kp . 
After identification of the connected clusters, the spin of 
cluster is flipped with probability 



pf P 



1 



1 + e -2a(C)"' ( 29 ) 

«(C)= ]T 9 p Pli ]T p(x)(3(x)p(x + ft). (30) 

j5—4>,a x£C,p 

Updating the j3 — <fi,a fields: Bonds between 
like-sign neighbors are kept with the probability 1 — 
e -2[Kf>-g P pfi e | —J, the spin of cluster is flipped with 
probability ^- 

In our simulation, the parameters are chosen as n p = 
0.3323, k p — 0.3897, n a = 0.3748, and g p ^ = g paa = 
0.02, the masses of (j> and a fields are measured through 
single particle propagators, the values are given by m^ ~ 
0.176 and m a ~ 0.240 respectively in lattice unit. The 
mass of resonance p is established from phase shifts, and 
the approximate value is given by m p ~ 0.57. 

In this work, we use T = 80 and various spatial exten- 
sions L between 15 and 50. For each set of lattice size 
and moving frame, we generated typically one million 
measurements. 



B. Particles spectrum 

As shown in the elastic scattering case in 1 + 1 dimen- 
sions |27j , one particle propagator can be constructed by 
operators 



^w^x;^)^ 1 "'". 



(31) 



0.19 — 
M f <L) 0.183- 
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FIG. 2: m^„ as function of L, they follow the function of 
M a (L) = m a + c a /L 1/2 e~ mc " L (red dashed curves). 



where q\. n = ^n, n = —L/2 + 1, • • • , L/2 and a — <f), a. 
The spectrum of single particle fields is extracted from 
exponential decay of the correlation functions 



(^a,n(X0J 



l (x )a„(0)) oc 



-K*q 



(32) 



The single particle's masses satisfy relation M a (L) — 
m a + c a L- x l 2 e- m <* L [21], see Fig{§ 

The two particles operators in the moving frame with 
total momentum of P = ^rd,d £ Z are constructed from 
single particle operators 

Of p , d) (x ) = p d (x ), (33) 

O( a , n )( x o) = a n (x )ad-n(xo). (34) 

The two particles correlation function matrices read 

Cg(xo) = ([Of* (x ) - S dfi Of*(x + 1)] Ofto)), (35) 

where short hand notation i,j denotes the different sets 
of (p,d) or (a,n). The disconnected contribution has 
to be subtracted in CM frame (d = 0). The spectral 
decomposition of the correlation function matrices has 
the form, 



i 

where v\ ' — (l\Of (0)|0) and I labels the energy eigen- 
state E) . The energy levels are determined by solving 



generalized eigenvalue problem [14 

C d (x )ti - X w) (x ,xo)C d (x ^i, 



(37) 



where \(d,l){xo,x ) 



-{x Q -x )E\' 



and Xq is a small 



FIG. 3: The energy spectra of a coupled-channel Ising 
model as function of (L,d): (Left) d — 0, (Middle) d — 1 
and (Right) d = 2. The red curves represent (I) the en- 
ergy spectra of a non-interacting pair of particles: E^ d ' = 
2~2i=± cosh -1 (coshm0, CT + 1 — cospi), where p± — ^n±± j-, 
n+ + n- — d, and (n±, d) £ Z. The masses of 4>, a are given 
by m$ ~ 0.176, me ~ 0.240; (II) the energy spectra of a sta- 
ble resonance: E^ d ' — cosh -1 (cosh m p + 1 — cos P), where 
m„ ~ 0.57. 



simulation, the size of the matrices varies according to the 
volume, the number of operators we are using is always 
two or three more than the number of energy eigenstates 
in the region 2m^, < y/s < Am$. The values of the energy 
levels are determined by fitting X^d,i)(xo,0) for $J xq ^ 
6—10 with the form 

A W) (so,0) = (l-Vi))^" 1,i ' ,IO +i W) e-'"w) I », 



where A 



(d,l)i m (d,l) 



and 



\d,i) 



are fitting parameters. 



This form allows a second exponential, however, we found 
that the value of m', d n is typically 2 — 3 times of the value 
of mid,l)i so that it decreases rapidly and the first expo- 
nential becomes dominant at around xq = 2. 

We show the measured two-particle energy spectra 
from our simulations in Fig|3]for various volumes and to- 
tal momenta of two particles system P = ^d, d = 0, 1, 2. 



reference time, in our analysis, xq is set to be zero. In our 



C. A coupled-channel _K-matrix model 

In order to extract the scattering amplitudes (phase 
shifts and inelasticity) from the discrete finite volume 
spectra of the Monte Carlo simulation, we consider a 
if-matrix model for a coupled-channel S-wave scattering 
system. 

In the scaling regime, the phase shifts of the Ising 
model (50. CT in 1 + 1 dimensions are shifted by a back- 



ground phase Si s i r , 



where the inverse of the t Res -matrix is given by 



<W - 5 6 *? - Sj, 



(38) 



where <5 Ros represent the normal phase shift in which a 
resonance may appear at value of S Rcs = \. Thus, the 
unitarized t-matrix may be defined by 



t aa = -C S + ^(s-4m^) 
t al 3 = -6{ S -Aml)tlf, 



Z/v™ 



where t— and £ Res -matrix are parametrized by phase 
shifts S<f, a , <5 R( ^ and inelasticity -q respectively The t- 



matrix is related to the/-matrix defined in Eq.(25l by 
equation t a g = — f= fad- The unitarity of the t-matrix 

is guaranteed by the unitarity of the £ Res -matrix. The 
Ising model suggestion is to parameterize a resonance 
coupling to both channels using a pole interfering with a 
polynomial in an S'-wave if-matrix, 



K a p{s) 



9a9P 

M 2 - . 



■<7 (0) 

la/3 



-7 (1) S- 
'a/3 b 



(39) 



J ap 



la/3 



Here I a (s) is the Chew-Mandclstam form [30 whose 
imaginary part above threshold (s > 4m„ ) is the phase- 
space, 



I a {s) = J a (0) - - / ds'xll 

T J Ami 



Am 2 ,, 



s' (s' — s)s' 



(41) 



We have opted to subtract the integral once, and it is 
convenient to choose I a (0) such that Re/ a (M 2 ) = so 
that we have an amplitude which for real s near M 2 is 
close to the Breit-Wigner form with mass M. 

Given an explicit model for the scattering amplitudes, 
we can solve Eq. ( 42 ) for the finite volume spectra in var- 



ious volumes and total momenta P — ^d, deZ. 




cot 



Pc/,L + nd 



Jv 



cot 



p a L + nd 



. p^L + ird \ I p, 

i + cot — — ) [ I + cot 



a L + Trd\ yv) 



f<P<P fa 
J-p J-p 



taa 
l-p J-p 



(42) 



Eq.(42| is derived from Eq.(24) by replacing ^k a with 
Pa (ct = 4't (t )i where p a is the relative momentum of 
two particles in a moving frame: p a = Pa ' 1 2 Pa ' 2 and 
P = Pa.i + Pa,2- To compensate for the ultraviolet cut- 
off effect from the finite lattice spacing, the dispersion 
relation cosh_E = coshrn + 1 — coshp [27 is used in this 



a function 



work, accordingly, in Eq. ( 42 ) , the relative momenta of 



two particles in a moving frame p a are solved by the 
equations 



P = Va,\ +Pa,2, 

E (d) ^ y^ cosh" 1 (cosh m a 



E 

i=l,2 



1 — cosp Q , 



The lattice dispersion relation and finite size effects are 
further discussed in Appendix [B] 



D. Data analysis 



With the if-matrix model described in Section |IVC 



we can perform a global fitting method proposed in [TJ 
to the spectra in Fig|3] For this purpose, we can minimize 



x 2 ({ 



E n (L,d) 



[E n (L,d)-E^(L,d;{a i })]' 
a(E n (L,d)) 2 



(43) 



within the space of if -matrix parameters, {a^} = 
{M,g$,g (r ,~f( n > . . .}, where E n denotes the energy levels 
from Monte Carlo simulation, and E^ ot are the solutions 
ofEq.@. 

Instead of establishing the resonance pole position in 
our toy model, the purpose of this work is to demonstrate 

(1) the methodology of extracting scattering amplitudes 
from data of coupled-channel Monte Carlo simulations, 

(2) predictability of scattering amplitudes extracted from 
a set of lattice data, and (3) the validity of our formalism 
while taking into account of the finite size effect presented 
in Appendix IB] Therefore, in this work, we choose to fit 
the spectra below 40 threshold for d — only, then for a 
consistency check, we compare our predicted spectra for 
d = 1, 2 to the spectra from the Monte Carlo simulation. 
We show the spectra of if-matrix model (red bands) in 
FigJ4] with the comparison of spectra from the Monte 
Carlo simulation (filled black circles, filled blue squares 
and filled green triangles) . The Zf-matrix we used in the 




FIG. 4: The finite volume energy spectra from /("-matrix 
model (red band) as function of (L,d), the spectra of K- 
matrix model are obtained by performing the fit on d = 
lattice data below 40 threshold only (black filled circles on 
the left). All the spectra in above three plots (d — 0, 1, 2) are 
presented in CM frame. The prediction of energy spectra from 
Tf-matrix model fit (red band) are also given for (Middle) d = 
1 and (Right) d — 2 compared to the Monte Carlo simulation 
data for d — 1 (blue filled square) and d = 2 (green filled 

triangle) respectively. -^ — = jy^-g = 0.93. 
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fitting has nine free parameters, the polynomial of 7'™- 1 
is taken up to 0(s 1 )(n = 0, 1). The value of parameters 
we find from fitting read 



M = 0.572(1),^ 
t2=0.3(1), 7 « 



0.064(4), g a = 

-0-7(3), 7^ 



0.060(4), 
= 0.11(3), 



7?J 



-0.3(i; 



1 1(7(7 



-0.6(2), 7^ = 1-5(5) 



0.48 0.5 0.52 0.54 0.56 0.58 0.6 0.62 0.64 0.66 0.68 0.7 



The extracted phase shifts S^,, 5 a and inelasticity r\ are 
shown in Fig(5j 

As demonstrated in the middle and the right plots in 
FigHJ our predicted energy spectra from if-matrix model 
(red bands) for d = 1, 2 agree with the spectra from the 
Monte Carlo simulation (filled blue squares and filled 
green triangles) within a reasonable precision. There- 
fore, we accomplished our goals, (1) we proved that our 
formalism gives the consistent result in different moving 
frames while taking into account of finite size effect, and 
(2) we have shown that the global fitting method is a 
valid and fairly reliable means for extracting scattering 
amplitudes from Monte Carlo simulation data. 



V. SUMMARY 

Based on the Lippmann-Schwingcr equation approach, 
in Section III] and |HI| we first derived a generalized 
Liischer's formula in 2D for two particles scattering in 
both the elastic and coupled-channel cases in moving 
frames. In Section |IV| we presented a 2D coupled- 
channel scattering lattice model. The model simulates 
a two-coupled-channel resonant scattering system, in 
which a resonance couples to both channels. Next, 
we performed Monte Carlo simulations on various fi- 
nite lattice sizes and in different moving frames. The 
discrete finite-volume spectra were extracted by fitting 
two-particle correlation functions. Finally, we used the 
2D generalized Liischer's formula to extract the scat- 
tering amplitudes for the coupled-channel system from 
the discrete finite-volume spectra. We have shown that 
the global fitting method can be used to reliably extract 
scattering amplitudes from Monte Carlo simulation data. 
The finite size effects on the solution of the generalized 
Liischer's formula were discussed in details in Section ITVl 
and Appendix[BJ We demonstrated that while taking into 
account of finite size effects, our formulae produce con- 
sistent results in different moving frames. 
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Appendix A: One dimensional infinite sum 

Let's consider the one dimensional infinite sum in 



Eq.(15l, 






i ik\x—x / — -fnL\ i-^nL 



(Al) 



In the region which we are interested in: \x\ > \x'\ and \x — x'\ < 7L, Eq.( Al I can be rewritten to 



n=£0 



e ik\x\ J2 Y v (x)Y v (x r )J^(kx') + J2 eikhnLle ^ nL J2 Y v (n)Y v (x)J v (kx) J^ Y v ,(n)Y v ,{x')J^{kx'). 

V=± n£Z V=± V'=± 

With the help of equations 

PL 



(A2) 



n#0 

J2e lkhnL] e^ nL Y + (n) 
nez 



cos ^ - e^ kL 
cos jkL — cos ^ 



n#0 

^e lfc|7 ™ L| e^" L r_(n) = i 



sin 



cos ~/kL — cos 



PL 



where the infinite sums are performed by using the property of polylogarithmic function IAq{x) — Xm=i x ™ 
thus, we find 



\ ^ ik\x — x — ~{nL\ innd 



Y^ Y v (x)Yp(x')J^,(kx') 



v=± 



J.k\x\ 



"fkL + 7rcT. 
I — icot — I J-p(kx) 



(A3) 



for |x| > \x'\ and \x — x'\ < yL. 

Appendix B: Lattice dispersion relation 

For the determination of the lattice spectrum at a pre- 
cise level, the finite size effect has to be taken into con- 
sideration by using the lattice dispersion relation 27J 



cosh^/i = cosh£ (d) - (1 - cosP), 



(Bl) 



r 



frame (fei + ki = 0) by Lorentz transformation relation 
p = jk. In the limit of vanishing lattice spacing, the 
Lorentz contraction factor 7 is given by 7 = —^. How- 
ever, due to the ultraviolet cut-off effect from finite lattice 
spacing, the definition of the Lorentz contraction factor 
7 = — — is inconsistent with lattice dispersion relation 
Eq.(Bl ). This inconsistency leads to the large discrepan- 



where E& and Js are the total energy of system in cies of P hase shifts and inelasiticy computed in different 
moving frames and the CM frame respectively, and the frame f ■ To resolve ^ thl ^£[ ^ len ^ e mayuse_the relation 
total momentum of system is given by P = jj_ d, d E Z. 
In the limit of vanishing lattice spacing, Eq.(Bl I reduces 
to the relativistic dispersion relation: E^ = \/s + P 2 . 



p = -fk to rewrite Eq. ( 19 1 , Eq. ( 24 ) and Eq. ( 26 1 to 



The relative momentum of two particles p 



, pL + nd 
cot dp + cot = 0, 



(B2) 



a moving frame (pi + pi = P) is related to the rela- 
tive momentum of two particles k = kl ~ k2 in the CM for single channel scattering, and 



1 p^L + nd 



1 p <T L + Trd\ ( . p^L + nd 

h I + COt = \ I + COt — 

f$ a 2 J \ 2 



i + cot 



PaL- 



nd\ (&)* 



f<P<P fo 

j-p Jt 



pa<7 
Pp Jv 



r 



Tfp{-1 



cos , 

,d \ 2 



**±!±L + 8++6?, 



cos 



-L + 8% - 6: 



■p u-p 



(B3) 



for coupled channel scattering, respectively. In Eq.(B2|, 
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FIG. 6: The energy spectra of single channel Ising model [23 
as function of (L, d) for (a) g = (Upper panel), (b) g = 0.02 
(Middle panel) and (c) g = 0.04 (Lower panel) respectively. In 
each panel, from left to right, each individual plot is related to 
d — (black filled circles) , 1 (red filled squares) and 2 (green 
filled triangles) respectively. The red curves represent (I) the 
energy spectra of a non-interacting pair of particles: E^ d ' = 
Ysi=± cosh -1 (coshm + 1 — cosp;), where p± — ^-n± ± f-, 
n+ + n- = d and (n±,d) £ Z; (II) the energy spectra of a 
stable resonance: E* ' — cosh -1 (coshm p + 1 — cosP), where 
m p ~ 0.5. 
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FIG. 7: The phase shifts of single channel Ising model [27] 
for (a) g = (Left), (b) g = 0.02 (Middle) and (c) g = 
0.04 (Right) respectively. The red line in left plot labels the 
theoretical expectation <5i S m g = — ? for non-interacting case 
(a) g = 0. The red filled circles, black filled circles, blue filled 
squares and green filled triangles represent the result in Fig. 6 
of [27], and our results for d = 0, 1, 2 respectively. 
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Eq. ( 42 ) and Eq. (|B3|) , the relative momentum of two par- 
ticles p = ^ 



is solved by equations 



P = Pi +P2, 
E {d) = Y^ cosh -1 (coshm + 1 - cosp,) . (B4) 



So that, rather than solving Eq.( 19 1, Eq.(24| andEq.(26) 



v^ ' 
we use Eq.(|B2|) and Eq.(|42|) with the solution of relative 



with the Lorentz contraction factor given by 7 



momentum of two particles given by Eq. ( B4 1 for single 



and coupled-channel scattering respectively in this work. 
As a simple demonstration how the above proposal 
works, let's consider a non-interacting two-particle sys- 
tem in 1 + 1 dimensions. Two particles in an arbi- 
trary moving frame, P = -^-d, have individual momenta 
p± = 75 n ± ^ J, respectively, where n + + n_ = d and 
{n±,d) € Z. So that we get a relation: P L + vd =n+7r+|, 
where p = p+ ~ p- is the relative momentum of two- 



particle system. Using Eq.(B2|, we conclude that the 



phase shift of two non-interacting particles is given by 



S = Rising = — f ■ This conclusion derived from Eq.(|B2|) 



holds in all the moving frames. 

As an more quantitative example, we generalize the 
single channel Ising model computed in CM frame in [27] 
to the moving frames, we compute the three sets of mod- 
els [37] (see Table 1 in [37J) for three different moving 
frames: d = 0,1,2. Three models are labeled by cou- 



10 



pling constants: (a) g = 0, (b) g = 0.02 and (c) g = 0.04, 
all the parameters are given by Table 1 in [27] . The mea- 
sured energy levels for three models (g = 0,0.02,0.04) 
and three moving frames (d = 0, 1, 2) are shown in Figpjj 
the extracted phase shifts by using Eq.(B2) along with 



the solution of relative momentum of two particles from 
Eq. ( B4 ) are shown in FiglTJ FigJ7] demonstrates the con- 
sistent calculation of phase shift from different moving 
frames by using Eq. (|B2|) along with Eq. ( B4 1 . 
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